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^ ; We analyze the possibility to distinguish between quintessence and phantom scalar field mod- 

els of dark energy using observations of luminosity distance moduli of SNe la, CMB anisotropics 
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and polarization, matter density perturbations and baryon acoustic oscillations. None of the 
present observations can decide between quintessence or phantom scalar field models at a statis- 
tically significant level: for each model a set of best-fit parameters exists, which matches all data 
with similar goodness of fit. We compare the relative differences of best-fit model predictions 
■ with observational uncertainties for each type of data and we show that the accuracy of SNe la 
luminosity distance data is far from the one necessary to distinguish these types of dark energy 
^ . models, while the CMB data (WMAP, SPT and Planck) are close to being able to distinguish 
them. Also a significant improvement of the large-scale structure data (e.g. Euclid or BigBOSS) 

in 
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will enable us to decide between quintessence and phantom dark energy. 



Introduction 

Scalar field models of dark energy are among the most promising and best elaborated ones to 
match observations of the accelerated expansion of the Universe. This is reflected in the numer- 
ous papers, reviews and text books on the subject (see for example [Ill2J[3j|IJ|5l[6l[71IEll9l [Hll EH] 
and citing therein). However, despite good agreement of theoretical predictions with modern 
observational data a number of problems remain unresolved and the prospects for their solu- 
tion are unclear. Among them is the question whether dark energy is quintessence, phantom, 
quintom, k-essence or of some other form. For each of these types of dark energy, cosmological 
best-fit parameters exist, for which the differences of maximum likelihoods are statistically in- 
significant for the same datasets. The identification of dark energy type is important in order 
to establish the properties of dark energy in our Universe, its physical nature, its origin and its 
couplings to other fundamental fields and particles. Without this knowledge, e.g., we cannot 
predict the future of our Universe. 

For example, in the paper [T2j [T5] we have shown that quintessence models of dark en- 
ergy with decreasing and increasing equation of state (EoS) parameter (freezing or thawing 
quintessence) cannot be distinguished by current observations, but could be so by future data. 
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In the recent papers [HI [15] it has been shown also that some datasets slightly prefer phantom 
models over quintessence, while in contrary other datasets prefer quintessence. However, in all 
the cases considered, the differences in the maximum likelihoods are not statistically significant. 

The goal of this paper is to estimate the possibility to distinguish between quintessence and 
phantom scalar field models using the current and future observational data. 



1 Scalar field dark energy and best-fit cosmological 
models 

We suppose that our Universe is spatially flat and filled with the non-relativistic particles (cold 
dark matter and baryons), relativistic particles (the thermal cosmic microwave background 
(CMB) and massless neutrino) and a minimally coupled scalar field with given Lagrangian. The 
scalar field is specified as follows: (i) its Lagrangian is canonical, L de = ±X - U(<f>) (X = <p 2 /2 
is the kinetic term and U(4>) is the potential), with "+" for a quintessence scalar field (QSF) 
and "-" for a phantom scalar field (PSF). The effective sound speed c 2 ^ = Spd e /Sp de (speed 
of propagation of the scalar field perturbations) is equal to the speed of light c for both cases; 
(ii) its equation of state is p de = Wd e c 2 Pde and the time derivative of the dark energy pressure is 
proportional to the time derivative of its energy density, p de = c 2 a c 2 f>d e - The coefficient c 2 , often 
called "squared adiabatic sound speed"0, is constant. Integrating this equation we obtain the 
generalized linear barotropic equation of state pd e = c 2 p de + C, where C is a constant, that is 
why we call such scalar field barotropic. 

We solve Einstein's equations for the background dynamics and the Einstein-Boltzmann 
system of linear perturbation equations in the synchronous gauge for the evolution of the 
perturbations. The background Universe is assumed to be the spatially flat homogeneous and 
isotropic with Friedmann- Robertson- Walker (FRW) metric, ds 2 = gijdx l dx^ = a 2 (r])(di] 2 — 
8 a pdx a dx^) ) where 77 is conformal time defined by cdt = a(r])drj and a(j]) is the scale factor, 
normalized to 1 today (below we set c = 1). We assume that dark energy and dark matter do 
not interact. The energy- momentum conservation for dark energy then determines Wd e and pd e 
as functions of the scale factor a, 

(l + cg)(l + Wo ) (0) (1 + w )a-^ + c 2 a -w 

Wd *- l + Wo _( WQ _ c 2) a 3(l+c3) Pde -Pde j + ^ ■ W 

Using the Friedmann equations we obtain the Hubble and deceleration parameters 



rr rr / o -4^0 3 T O f( \ 1 2^0"* + O^ 3 + (1 + 3 W<fe )£W(a) 

H = H ^/Q r a 4 + Q m a 3 + Q d ef{a), q = — _ 3 

2 \l r a 4 + iz m a 6 + life j [a) 



(2) 

where f(a) = p de / p^ , p d °] is current value of dark energy density, w is EoS parameter Wd e 
today, H is current value of Hubble parameter (Hubble constant) and Q r = Pr / Ptot-i = 
Pm / Ptoti ^de = Pfe / Ptot are ^ ne dimensionless density parameters of the relativistic, non- 
relativistic and dark energy components correspondingly. The first Firedmann equation requires 
Q r + fl m + Qde — 1- The matter density parameter is the sum of cold dark matter and baryons, 

\L m = licdm + i'b, and p tot = p r + pm + p de ■ 



1 It corresponds only formally to the adiabatic sound speed in thermodynamics. 
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Since the light neutrinos have become non-relativistic by today, the density parameter 
of the present relativistic component is given only by the photon density which is accurately 
determined by the value of the CMB temperature, fi r = fi 7 = lQirGa SB To /3H$ = 2.49 -lO" 5 /*- 2 
and can be neglected in the matter and dark energy dominated epochs. The density parameters 
of the other components are less well known and depend somewhat on the model of dark 
energ)jE The scalar field affects the expansion of the Universe, it causes accelerated expansion 
when |(1 + ?>Wde)^def{,o)\ > f2 m a~ 3 as follows from eqs. (j2j). 

Using expressions (CQ) and (j2|) we can compute the "luminosity distance - redshift" or "an- 
gular diameter distance - redshift" relations to determine all the above-mentioned parameters 
by comparison with corresponding observational data on standard candles (supernovae type 
la, 7-ray bursts or other) and standard rulers (positions of the CMB acoustic peaks, baryon 
acoustic oscillations, X-ray gas in clusters or other). 

We assume the standard paradigm of large scale structure formation: (i) it is formed from 
stochastic, adiabatic, Gaussian scalar perturbations generated in the early Universe, (ii) the 
initial power spectrum of radiation and matter density perturbations is power-law, Pi{k) = 
A s k ns , where A s and n s are the amplitude and spectral index (k is wave number). The scalar 
field cannot be perfectly smooth, it is perturbed by gravitational influence of matter-radiation 
inhomogeneities or has its own initial fluctuations, generated in the early Universe. 

The system of linear differential equations for the evolution of quintessence and phantom 
scalar field perturbations and their numerical solutions are analyzed in our previous papers 
[HI EH E]. The main conclusions are as follows: (i) the amplitude of scalar field density 
perturbations at any epoch depends strongly on parameters of barotropic scalar field fide, wq, 
c a and c 2 ; (ii) although the density perturbations of dark energy at the current epoch are 
significantly smaller than matter density perturbations, they leave noticeable imprints in the 
matter power spectrum, which can be used to constrain the scalar field parameters. 

The linear power spectrum of each component can be computed by numerical integration 
of the Einstein-Boltzmann equations [HI [191 1201 [21] using publicly available codes such as 
CMBFAST [221 [23], CMBEasy [21], CAMB [23 EE] or CLASS [23 EEl EE] , 

Piin(k) = Pi(k)T 2 (k] Q r , Q b , Q cdm , Q de ) . 

T(k) is the transfer function, for the cosmological model with the given parameters. In all the 
computations presented here and in recent papers, we use CAMB with corresponding modi- 
fications for the quintessence/phantom barotropic scalar field as dark energy. Comparison of 
computed matter density power spectra with the ones obtained from galaxy surveys constrains 
the parameters of the models discussed here. 

The same Einstein-Boltzmann equations and codes also determine the angular power spec- 
tra of CMB temperature anisotropics, Cj T , and polarization, Cf E , as well as their correlations, 
Cf T , which can be compared with WMAP data to constrain the cosmological and DE param- 
eters mentioned above. The calculation of CMB anisotropics and polarization requires also the 
knowledge of the reionization history of the Universe, which depends on complicated non-linear 
effects of structure formation. This is parameterized by the value of optical depth from current 
epoch to decoupling (at redshift Zd ec ) caused by Thomson scattering. It is denoted by Td ec and 
is also fitted by the data. As it is supported by numerous papers (see for example reviews 
[DElElllElElElElISlIinKIIJand citing therein), the CMB data are the most important ones 
in the determination of cosmological parameters. The parameter estimation from the data is 
performed using publicly available Markov chain Monte Carlo (MCMC) codes [30[ |3T]. 

2 Hcre and below h = i7 /100km/sMpc 



3 



Table 1: The best-fit values p_f] and approximate 2a marginalized confidence ranges of cos- 
mological model parameters for QSF (q) and PSF (p) determined by the MCMC technique 
using 2 observational datasets: WMAP7 + HST + BBN + BAO + SN SDSS SALT2 ( qi , Pl ) 
and WMAP7 + HST + BBN + BAO + SN SDSS MLCS2k2 (q 2 , p 2 ). We denote the rescaled 
energy density of the component X by ux = flxh 2 . 
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The cosmological model with quintessence/phantom barotropic scalar field with canonical 
Lagrangian has q free parameters f2&, fl c dm, Qde, wo, c 2 , H , A s , n s , Td ec in the general case and 
8 when we restrict ourselves to spatially flat cosmological models, since fi& + ^cdm + ^de = 1- 
To constrain these parameters we used the following datasets: 

1. CMB temperature fluctuations and polarization angular power spectra from the 7-year 
WMAP observations (hereafter WMAP7) [3_l __B _3] • 

2. Baryon acoustic oscillations in the space distribution of galaxies from SDSS DR7 (here- 
after BAO) [35]; 

3. Hubble constant measurements from HST (hereafter HST) [36] . 

4. Big Bang Nucleosynthesis prior on baryon abundance (hereafter BBN) [371 138] . 

5. Supernovae la luminosity distances from SDSS compilation (hereafter SN SDSS) [39|, 
determined using SALT2 method of light curve fitting PS] (hereafter SN SDSS SALT2) 
and MLCS2k2 [41] one (hereafter SN SDSS MLCS2k2). 

In our previous paper p_t] we have performed the Markov chain Monte Carlo (MCMC) analysis 
for two combined datasets: WMAP7 + HST + BBN + BAO + SN SDSS SALT 2 (dataset 
1) and WMAP7 + HST + BBN + BAO + SN SDSS MLCS2k2 (dataset 2) to determine the 
best-fit values and confidence limits of the model parameters for QSF/PSF. The results are 
presented in Table [U The upper/lower limits are approximate marginalized 2a ranges from the 
ID posteriors. In the last row the negative of the logarithms of maximum of likelihood functions, 
X 2 = —log(L max ), are shown for comparison. The PSF+CDM model with best-fit parameters 
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p l determined on the base of dataset 1 has a slightly lower x 2 than the QSF+CDM model 
with best-fit parameters q : . In the case of dataset 2 the result is opposite: the QSF+CDM 
model with best-fit parameters q 2 matches the dataset 2 better than the PSF+CDM model 
with best-fit parameters p 2 . 

In Ref. [14] we have also determined the best-fit parameters of ACDM from the datasets 

1 and 2 and we have found that in both cases the x 2 f° r ACDM is between the correspond- 
ing values for QSF+CDM and PSF+CDM. This trend suggests that the best-fit parameters 
of QSF+CDM and PSF+CDM (except for c 2 ) are determined reliably. But in both cases the 
differences of the x 2 's are not statistically significant. Hence, we have a degeneracy between the 
QSF and PSF models of dark energy and we shall elucidate here how strong this degeneracy 
is and whether the accuracy of current or forthcoming data is sufficient for lift it. Observa- 
tional constraints for another class of scalar fields with barotropic equation of state without 
peculiarities in the past are discussed in Appendix [A] 

The results for the cosmo logical parameters, especially H , Q de , Wd e and c 2 a , presented in 
Table HJ also indicate a tension between the two fitters SALT2 and MLCS2k2 applied to the 
same SNe la. This has already been highlighted and analyzed in Refs. [39l 14*2] . but up to now 
we have no decisive arguments for favor of one of the two lightcurve fitters. 

2 Distinguishing quintessence and phantom field models 
of dark energy 

The QSF+CDM model with best-fit parameters and PSF+CDM model with best-fit pa- 
rameters Pj all provide a good fit to the data from SNe la distance moduli (left panel of Fig. 
[1]), BAOs (left panel of Fig. [3]), matter power spectrum (left panel of Fig. [4j, CMB temper- 
ature fluctuations, polarization and cross-correlations polarization (left panel of Fig. [5]). Let 
us analyze the possibility to distinguish between QSF+CDM and PSF+CDM by each type of 
data. 

The lines corresponding to different DE models in the left panel of Fig. [T] look perfectly 
superimposed. The maximal values of relative differences of SNe la distance moduli (|/i(qj) — 
/i(pi)|/yu(qi)) for QSF+CDM and PSF+CDM models with best-fit parameters determined using 
the same fitters are less than 0.1% (right panel of FigJT]). Comparing this to the observational 
uncertainties plotted as dots and triangles, we conclude that these data are very far from 
distinguishing between the proposed types of scalar field dark energy. 

Fortunately, other characteristics of the expansion dynamics of the Universe, based on 
measurements of the first and second time derivatives of the Hubble parameter H(z), are 
significantly more sensitive to the current value and time dependence of the EoS parameter. 
The redshift dependence of dimensionless parameters describing the expansion dynamics of the 
Universe, such as the rate of expansion H/Hq, the deceleration parameter q = —H/(aH 2 ) — 1 
and the statefinder parameters [43] 



are shown in Fig. [2] for the models with best-fit parameters qi, q 2 , Pi and p 2 . In our dark 
3 /J,(z) = m — M = 51ogc?L + 25 is the integral of 1/H over redshift, since cLl{z) — (1 + z)c J* dz' / H(z'); 



r = H/(a 2 H 3 ) + 2H/{aH 2 ) + 1, s = (r - l)/3(g - 1/2) 



(3) 



H(z) can be deduced from observations by differentiating of fi(z) with respect to z, since H(z) 
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Figure 1: Left panel: the distance modulus fi = m — M as function of the redshift z for SNe la 
in the models with best-fit parameters qi, q 2 , Pi and p 2 (superimposed lines) are compared to 
the data from SDSS SNe la (symbols). Right panel: the relative differences of distance moduli 
in QSF+CDM and PSF+CDM with best-fit parameters determined from the same datasets 
(lines) are compared with observational uncertainties for SDSS SNe la (symbols). 
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Figure 2: The redshift dependence of dimensionless parameters describing the expansion dy- 
namics of the Universe for the models with best-fit parameters qi, q 2 , pi and p 2 . The rate of 
expansion H/Hq (left panel, top), the deceleration parameter q (left panel, bottom) and the 
statefinder parameters r and s (right panel) are shown. The dashed ranges are their observa- 
tional uncertainties estimated by variation of the expressions d2J) for H/H , q, (J3J) for r and s 
over Q de , w , c 2 a for 2a marginalized uncertainties of all models presented in Table [TJ 
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Figure 3: Left panel: the relative BAO distance measure R r D = r s (zd rag )/ 'Dv(z) for the mod- 
els with best-fit parameters qi, q 2 , pi and p 2 (lines) is compared with data from the SDSS 
DR7 [32] and WiggleZ [44] surveys (symbols). Right panel: the relative differences of the BAO 
distance measure \AR rD \/ R rD in the models with best-fit parameters qj and p { . Symbols show 
observational la errors of the data in the left panel. 

energy models ([1]) the state finder parameters reduce to 

r = 1 + 4.5(1 + w de )c 2 a tt de (a), s = (1 + w de )c 2 a /w de , (4) 

where fl de {a) = 87iGp de (a)/3H 2 . The differences between H/ H for q^ and p i; as well as between 
g(q.j) and g(pj, s(qj and s(p i ), r(q i ) and r(pj, at different z are significantly larger than those 
of fi(z) or (1l(z). Unfortunately, current data on cIl{z) from SN la measurements are too poor 
to determine these quantities. The dashed ranges in Fig. [2] show the dispersion obtained by 
the variation of expressions ([2} for H/H and q, (jlj) for r and s over Q de , Wq, c a within the 2a 
uncertainties of all models presented in Table [TJ Maybe future high-precision measurements 
of distances to a significantly larger number of SN la or extremely well-calibrated 7-ray bursts 
will provide a possibility to distinguish the QSF and PSF models of DE. But clearly, present 
SNe la data is very far from this goal. 

The relative BAO distance measure Rtd(z) = r s (z drag )/Dv(z) (r s (z drag ) is the sound horizon 
at the epoch where photon drag stops due to recombination and Dy{z) is the BAO dilation 
scale) extracted from SDSS DR7 galaxy redshift survey [35] and from the WiggleZ survey [44J 
is matched well by the QSF+CDM and PSF+CDM models with best-fit parameters qi, q 2 , 
Pi and p 2 (see Fig. [31 left panel). The relative differences of R t d(z) for models with QSF and 
PSF are up to 2% for z < 0.4, while observational errors are about 3% at 0.2 < z < 0.35 and 
more than 4.5% at 0.44 < z < 0.73 (right panel of Fig. [3]). This means that an improvement 
by a factor 2 to 3 in the determination of R r r>(z) in future large scale structure surveys like 
BigBOSS [4S] or Euclid [4B] can distinguish between these types of dark energy. The power 
spectrum of matter density perturbations extracted from luminous red galaxies sample from 
SDSS DR7 catalogue by Reid et al. (2010) [47] has not been used for the determination of best- 
fit parameters of QSF+CDM and PSF+CDM, however the computed power spectra match it 
well too. The experimental errors of its determination are still too large (8 — 12%) to distinguish 
between different scalar field models of dark energy. In addition, at the small scales (k > 0.1 
hMpc -1 ) there are uncertainties in the computation of the power spectrum associated with the 
nonlinear evolution of perturbations. The resolution of this problem would require combined 
N-body + scalar field simulations. 
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Figure 4: Left panel: the power spectrum of matter density perturbations in the cosmological 
models with best-fit parameters qi, q 2 , pi and p 2 . Dots show the observational SDSS LRG 
DR7 power spectrum [27]. Right panel: the relative differences of matter density power spectra 
\AP(k)\/P(k) in the models with best-fit parameters and p^. The dots show the observational 
uncertainties (la) of SDSS LRG DR7 data [27] . 



Above we have discussed the importance of CMB data for the determination of cosmolog- 
ical parameters and, in particular, for dark energy parameters. The key cosmological data at 
present are the WMAP all-sky maps, which contain information about the primordial fluc- 
tuation amplitude and spectral index, the positions and amplitudes of the acoustic peaks, as 
well as the amplitude of large scale matter density perturbations at late time imprinted on the 
integrated Sachs- Wolfe effect. In Fig. [5] (top left panel) the power spectra of temperature fluc- 
tuations extracted from the 7- year [32j [33j [32] and 9- year [28], |29] WMAP all-sky measurements 
are shown. Its accuracy is best (minimal errors ~ 1.5 — 4%) in the range of the first and second 
acoustic peaks {£ ~ 200 — 600). This allows an accurate determination of the main cosmological 
parameters. The best accuracy of dark energy parameters is achieved when CMB data are used 
together with SNe la data and relative BAO distance measures or directly the matter power 
spectrum. The power spectra £(£ + l)Cf T /2n for the QSF+CDM and PSF+CDM models with 
best-fit parameters qi, q 2 , pi and p 2 provide good fits to the WMAP7 and WMAP9 data, as 
well as to other recent CMB experiments, ACT J5U] and SPT [SI], which are sensitive to higher 
values of £. The relative differences between the power spectra in QSF+CDM and PSF+CDM 
models, shown in the top right panel of Fig. [5j do not exceed the relative observational uncer- 
tainties of £(£ + 1) Cj T /2tt for ACT. The relative difference between the spectra in models with 
parameters q 2 and p 2 is comaparable to the relative uncertainties of WMAP7 and especially 
WMAP9 data at £ ~ 300 — 800 and to the relative uncertainties of SPT data at higher £. 

We also present the forecasted uncertainties of the Planck data (based on the sensitivities 
for 7 frequency channels which can measure the polarization [52J, f s k y is taken to be 0.65). For 
I < 500, the WMAP errors are already nearly cosmic variance limited and the Planck data will 
not bring a significant improvement. At high £, however, the forecasted Planck uncertainties 
are smaller (for q 2 and p 2 ) or comparable (for q x and pi) to the relative differences of models. 
This is in good agreement with our previous conclusion [T2~] [T3] that the accuracy of the power 
spectra expected from Planck will significantly narrow the allowed range of parameter values 
for scalar field models of dark energy. 

Additional constraints on cosmological parameters are obtained when CMB polarization 
data are included. To illustrate the agreement between theory and observations, the temperature- 
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Figure 5: Left panel: the temperature (TT), temperature-polarization (TE) and polarization 
(EE) CMB angular power spectra, £(£ + l)Cf T /2ir, (£ + l)Cf E /2n and t{t + l)Cf E /2n, for 
the cosmological models with best-fit parameters qi, q2, pi and P2 (superimposed lines) are 
compared to currently available data (symbols). Right panel: the relative differences of tem- 
perature \ACj T \/Cj T and polarization \ACf E \/Cf E power spectra and absolute differences 
for the temperature-polarization correlation, (i + l)\ACf E \/27r, in the models with best-fit 
parameters q^ and (Table [1]). The symbols show la uncertainties of currently available and 
forecasted future data. 

polarization [I + l)Cj E /2n and polarization l{t + l)Cf E /2ti CMB power spectra for the cos- 
mological models with best-fit parameters qi, q2, pi and P2 are presented in the middle and 
bottom left panels of Fig. [5j All lines are superimposed for i > 20 with sub-percent accuracy 
for TE and a few percents for EE, while the errors of the WMAP7 and WM AP9 power spectra 
at these E's are larger than ~ 3 — 6% for TE and ~ 20 — 40% for EE, as it is shown in the 
middle and bottom right panels of Fig. |5j In the bottom right panel of Fig. [5] we present 
also the relative uncertainties from the QUIET experiment [53J (for both maximal likelihood 
and pseudo-C^s). The minimal errors of the QUIET EE power spectra at high I are at least 
~ 10 — 20%. The forecasted relative uncertainties for Planck are comparable to the differences 
between the models at high £. All this means that the data on polarization of CMB are unsufli- 
cient for the problem of distinguishing between quintessence and phantom fields at the present 
level of accuracy of observational data, but will possibly become sufficient in near future. 
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3 Conclusion 



None of the used observational data at the current level of accuracy can prefer quintessence 
cosmology (QSF + CDM) or phantom cosmology (PSF + CDM) at a statistically significant 
level. In the framework of each model a set of best-fit parameters exists. The best-fit model 
matches well each type of data and all together with similar goodness of fit. At present, the 
data on CMB temperature fluctuations from WMAP9 and SPT and the expected Planck data 
look most promising for the purpose of distinguishing between these models of dark energy. 
In the future, increasing the accuracy of CMB fluctuation measurements jointly with high 
precision matter density data will probably give the possibility to establish the dynamical 
properties of dark energy and, maybe, its nature. On the other hand, SNe la luminosity 
distance measurements do not seem very promising. 
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A Observational constraints on scalar field models with 

Wo > — 1, c 2 a < — 1 

Another subclass of scalar field models of dark energy with barotropic equation of state without 
peculiarities in the past are fields with wq > —\, c 2 a < —\. In the past the EoS parameter of 
such fields evolved from —1 at early epoch to wq today. The properties of these models will be 
studied in detail in a separate paper, here we only present observational constraints on them 
from the datasets 1 and 2, determined by the MCMC technique in the same way as in [14]. 

In Fig. [6] we see that while for wq the posteriors and mean the likelihoods are close and 
have the shape of Gaussian and half-Gaussian, for c 2 the shapes of the likelihoods differ signif- 
icantly from each other and from Gaussian. This means that also for such fields the range of 
acceptable values of c 2 is generally unconstrained (it is constrained only in the case q2). The 
ranges of acceptable values of all other parameters of the models studied here including wq 
are well constrained. In Table [2] we present the best-fit values and approximate 2a marginal- 
ized confidence ranges of cosmological parameters in models with scalar fields with wq > — 1, 
c 2 < —1. The value of x 2 f° r this model in case of dataset 1 is larger than the corresponding 
value for both quintessence and phantom dark energy. On the other hand, for the dataset 2 the 
value of x 2 is smaller than for the quintessence and phantom models. Thus the scalar field with 
w > —1, c 2 a < —1 is preferred by the dataset 2, however, again the difference is not statistically 
significant. 

Comparing the relative differences of distance moduli, BAO relative distance measures, mat- 
ter density perturbations and CMB power spectra between such models and either quintessence 
or phantom we see that scalar fields with w Q > —1, c 2 < —1 could only become distinguished 
from them in near future by Planck data. The probability of distinguishing them from phantom 
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Figure 6: One-dimensional marginalized posteriors (solid lines) and mean likelihoods (dotted 
lines) for wo (top panels) and c\ (middle panels) in models with u> > — 1, c 2 a < — 1. Left: 
WMAP7 + HST + BBN + BAO + SN SDSS MLCS2k2. Rright: WMAP7 + HST + BBN 
+ BAO + SN SDSS SALT2. Bottom: the corresponding two-dimensional mean likelihood 
distributions in the plane c\ — Wq. Solid lines show the la and 2a confidence contours. 

fields is slightly higher than from quintessence ones. At current level of accuracy the models 
are indistinguishable by all data. 
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